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Abstract 

A previous approach with Fermi-Dirac distributions for fermion partons is here 
improved to comply with the expected low x behaviour of structure functions. We 
are so able to get a fair description of the unpolarized and polarized structure 
functions of the nucleons as well as of neutrino data. We cannot reach definite con- 
clusions, but confirm our suspicion of a relationship between the defects in Gottfried 
and spin sum rules. 
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1 Introduction 



The experimental data on the unpolarized and polarized structure functions of the nucle- 
ons suggest a role of Pauli principle in relating the shapes and the first momenta of the 
distributions of the various quark parton species, for which Fermi-Dirac distributions in 
the X variable have been proposed. Here, to comply with the expected behaviour for the 
structure functions in the limit x — > 0, we add a liquid nnpolahzed component dominating 
the very low x region and not contributing to quark parton model sum rules (QPMSR). 
Section 2 deals with the motivations for our description of parton distributions and the 
comparison with available data. The results found and the implications for QPMSR are 
discussed in section 3. Finally we give our conclusions. 

2 Pauli exclusion principle and the parton distribu- 
tions 

The behaviour at high x of F2{x)/F2{x) [U, known since a long time, and the more recent 
polarization experiments [@, 0, which show that at high x the partons have spin parallel 
to the one of the proton, imply that u^x) is the dominating parton distribution in the 
proton at high x. 

Indeed at = 0, the axial couplings of the baryon octet are fairly described in terms of 
the valence quarks 

''li = i + F uL,^l-F , (1) 

C.i±£^ . (2) 

With the actual values [H, p 

F = 0.464 ±0.009 D = 0.793 ± 0.009 , (3) 

ul^i is larger that the others 



u 



val 



In a previous work we assumed that the parton distributions at a given large depend 
on their abundance at = 

p{x) = T{x,pyai) , (5) 



with T an increasing function of 'p^ai and with a broader shape for higher values of "p^ai- 
This assumption and the observation that with F = 1/2 and D = 3/4 (very near to 
the values quoted in Eq.(^)) one has u\j^i = 1/2 just on the center of the narrow range 
(4a/,4aJ = (3/8,5/8) lead to take 



1 r 



1 



u\x) = - S{x)^d^[x) =-d[x) , (6) 
which implies 

AM(a;) = — = u(a;) — (i(a;) . (7) 

Eq.(^ connects the contribution of IS.u[x) to the polarized structure function of the 
proton g\{x\ with the contributions of u and d to the difference between the proton and 
the neutron unpolarized structure functions F2{x) 



xg\{x] 



Since a smaller negative contribution to g\[x) is expected from I\d[x) for the twofold 
reason that = (l/4)e^ and IS^d^ai ~ —{Ij^/Suyai-, Eq-® should hold with a good 
approximation for the total structure functions in the region dominated by the valence 
quarks: this is just the case for x > 0.2 [0. 

By integrating Eq.(|^ one relates Am, which contributes to the spin sum rules, to u — d 
contributing to the Gottfried sum rule 



/ -[F^{x)-F,^{x)] dx = -{u + u-d-d) = - , (9) 
Jo X 3 3 



where the last equality follows if the sea is SU{2)i invariant {d = u). Indeed NMC 
experiment gives for the l.h.s. of Eq.(|^) 

Jg = 0.235 ± 0.026 , (10) 

implying 

J - M = 0.15 ±0.04 M - (i = 0.85 ± 0.04 . (11) 



Many years ago Field and Feynman suggested |]T0|] that Pauli principle disfavours the 
production of uu pairs in the proton with respect to dd, since it contains two valence u 
quarks and only one d. The correlation shape-abundance for the parton distributions is 
also the one suggested by the Pauli principle: the most abundant parton is the one 



dominating at high x and the assumption that and d have about the same shape seems 
confirmed by the experiment. 

The role of the Pauli principle has suggested to assume Fermi-Dirac distributions in 
the variable x for the quark partons [0] 



p{x) = f{x) 



fx — x(p)\ 
exp I + 1 



X 



(12) 



where f{x) is a weight function, x plays the role of the temperature and x{p) is the 
thermodynamical potential of the parton p, identified by its flavour and spin direction. 
Consistently for the gluons, neglecting their polarization, one assumed 



G{x) = ^fix) 



x-xg , ^ 
exp I I - 1 



(13) 



where the factor 16/3 is just the product of 2 {Sz{G) = ±1) times 8/3, the ratio of the 
colour degeneracies for gluons and quarks. To reduce the number of parameters, the 
distributions of d, s and of their antiparticles have been given in terms of the ones for u 
and u 

u\x) 



d{x) 

Ad{x) 
d\x) 
s{x) 
As(x) 



1-F ' 

—k f{x) exp I 

d^(x) = u^{x) 



X — x(u^) 



X 



exp 



X — X[U 



X 



1 



1 



six) 



As(x) 



u{x) + d{x) 







(14) 

(15) 
(16) 
(17) 
(18) 



with k fixed by the condition Ad = Ad^ai = F—D. For the weight function one considered 
the simple form f{x) = Ax°'. In terms of seven parameters, x, the x for u and u, A and 
a, one obtains a nice description of the unpolarized and polarized structure functions 
for the nucleons, but it is not possible to reproduce the the fast increasing |11| at low 
X of q{x) = u{x) + d{x) + s{x), confirmed at very high from the behaviour of -Ffl^) 
measured at Hi |1T2| . 

Indeed, we know that the form given in Eq. ([T2|) , with different values of the x for the 
different parton species, is not suitable in the limit x — > for the most divergent part 
expected on general grounds to be equal for the different partons, at least in the limit 
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of flavour symmetry. This most divergent part should not contribute to QPMSR as the 
ones given by Gottfried and Bjorken with 1 = 1 quantum numbers exchanged. It is 



therefore needed, to reproduce the data and to get information on the status of QPMSR 
within this approach, to disentangle the most divergent part for x ^ of the parton 
distributions (which does not contribute to QPMSR) from the remaining part, just given 
by Eq.dg). 

To this extent we add a liquid unpolarized component, giving to the light partons, u, 
d and their antiparticles, the same contribution = Al — x)^^, and to s and s, 

as in Eq.(0), fL{x)/2. To be not influenced by theoretical prejudices we consider as free 
parameters the x for u^^^\ d^^^\ u and d. Finally we introduce a new parameter in /(x) 

f{x) = Ax^'il-xf . (19) 

We try also to describe the structure function 

F^{x) = u{x) + d{x) + s{x) — u{x) — d{x) — s{x) , (20) 

measured with great precision in deep inelastic reactions induced by (anti)neutrinos [ p^ . 
According to Eq. (P^D , we expect s{x) — s{x) = and thus -Ff (x) — F^i^x) and -F3(x) to 
depend on the differences u{x) — d{x) and d{x) — u{x). Note that we cannot impose the 
conditions 

u-u = 2 , (21) 
d-d = 1 , (22) 

since they would imply the Gross and Llewellyn-Smith sum rule [15 

^F3(x)dx = 3 , (23) 
Jo 



which as well-known experimentally shows a defect. The l.h.s. of Eq.(p3D is in fact mea- 
sured to be 2.50 ± 0.018 (stat.) ± 0.078 (syst.) [|T^, defect commonly explained in terms 
of QCD corrections |T^. 

In Figures 1-6 we compare our predictions for -Ff*!^) ~ -^2* (2^)5 ^^si 

) and 

X(7"(x), which do not receive contributions from the /zguzc? component, and for F^(x)/F2'(x) 
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and xq{x) with the experiments. We restrict the x's to be < 1, since the factor (1 — x)'^ 
in f{x) makes the dependence of the distributions on x > 1 very smooth. 

Our initial goal was to introduce spin-dependent x's also for the g's and to test the 
relationship 

Au{x) = u{x) - d{x) , (24) 
assumed in previous works (0] and Eqs.(0) and (p^) would imply the equality 



xgf{x] 



Unfortunately, we found practically the same with negative and positive values for 
Au{x) and/or Ad{x), and realized that, with f{x) to be found from the data, we are 
unable to disentangle the contributions of Aq{x) and Aq{x) to the polarized structure 
functions. Thus, our choice Au{x) = Ad{x) = is neither motivated by data, nor by 
theoretical prejudices, but simply from our present inability to get information even on 
their signs and to settle the important issue, relevant also for the validity of the Bjorken 
sum rule, whether Eq.(^4|) is satisfied. Eq. (|25|) , with only An, Ad, Au and Ad contributing 
to the polarized structure functions of the nucleons, would imply 

5 



X 



mx)-F-{x)] . (26) 



In Table 1 we report the parameters found here by means of the MINUIT fitting code, as 
well as the ones of the previous fit (without liquid) and of a one by Bourrely and Soffer 
found on similar principles, but with several different assumptions from ours. 

Indeed, the data on unpolarized nucleons structure functions are at = 4 GeV"^ , 
the neutrino data at = 3 GeV^ |]T3[, and q measures are performed at = 3 GeV^ 
and 5 GeV^ |TT]] and differ at small x, while our curve is intermediate between the two 
sets of data. The data on gi{x) are at = 2 GeV^ [|l8l, whereas g\{x) is measured at 

= 10 GeV^ by SMC [| and at = 3 QeV^ by E143 §; despite some narrowing of the 
distribution at higher showing up in the data, the values of Ip are in good agreement. 
This fact and the expected dependance |T^, smaller than the actual errors on the 
polarized structure functions, gives us confidence that our analysis is shghtly affected by 
our neglecting the dependance of the distributions. 
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The parton distributions so found are described in Figure 7. The total momentum 
carried by q and q is 53%. In order that gluons carry out the remaining part of the mo- 
mentum xg is fixed to be —1/15. The gluon distribution is compared with the information 
found on them in CDHSW P|, SLAC+BCDMS m and in NMC pTf experiments at 



Q = 20 GeV in Figure 8. The agreement is fair for x > .1, while the fast increase at 



small X, confirmed also from the data at very small x at Hera |]22|, confirms that a liquid 
component is needed also for gluons. The excess at high x of our curve with respect to ex- 
periment may be, at least in part, explained by the expected narrowing of the distribution 
from = 4 GeV^, where we fit the unpolarized distributions, to = 20 Gey^Q 

3 Discussion of the results 

The inclusion of the liquid term and the extension of our fit to the precise experimental 
results on neutrinos has brought to substantial changes in the parameters with respect to 
the previous work 0. 

The low X behaviour of f{x) become smoother 2;--203±.oi3 instead of but this 

is easily understood since the previous behaviour was a compromise between the smooth 
gas component and the rapidly changing liquid one to reproduce the behaviour of q{x). 
The liquid component, relevant only at small x, carries only .6% of parton momentum 
and its behaviour ~ x~^'^^, similar to the result found in |2^, is less singular than the 
one, suggested in the framework of the multipherial approach to deep-inelastic scattering, 
proportional to ~ x^^'^ [^. The parameter x{u^) took the highest value allowed by us 
(1.), since the factor in f{x), (1 — x)^'^^, is taking care to decrease u\x) at high x. The 
temperature x is larger than the previous one and the one found by Bourrely and Soffer 
1% . Instead x{u^) is slightly smaller than the previous determination and about half 



the value found in |[T^, where f{x) is different for and uK 

The ratio r = u^{x)/d{x) varies in the narrow range (.546, .564) in fair agreement with 

the constant value 1 — F = .536 ± .009 assumed in and slightly larger than the value 

* Indeed the gluon distributions are obtained from the dependance of the distributions according 
to the LAP equations |2^. In this respect it is worth noticing that the parameter Aqcd found from the 
corrections to the scaling is slightly smaller than the one found from different sources Q . This qualitatively 
supports the idea that the evolution equations may be modified as a consequence of quantum statistical 
effects which would favour harder quarks and softer gluons, giving rise to a slower softening of quark 
distributions with increasing 
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1/2 taken in § and 0. 

The central value found for the first moment of Ugasix), .03, is smaller than dgas{x)/2, 
.08, while Eq.(|2^) implies u{x) > d{x)/2. However, the large upper error on Ugas and the 
uncertainty in disentangling the gas and liquid contributions for the g's do not allow to 
reach a definite conclusion about the validity of Eq. (p4D . 

Indeed our distributions are very well consistent with Eq. (PBD , as it is shown in Figure 9, 
where our predictions for the two sides of Eq. (^6|) are compared. 

We have been suggested by Prof. Jacques Soffer to compare the parton distributions 
found here with the measured asymmetry for Drell-Yan production of muons at y = in 
pp and pn reactions 

^ d(Tpp/ dy - daprj dy 
dapp/dy + dapn/dy 

which at rapidity y = reads 

^ jXsix) - l)(4A(x) - 1) + (A(x) - l)(4A,(x) - 1) 
"""^ (A,(x) + l)(4A(x) + l) + (A(x) + l)(4A,(x) + l) ' ^ ^ 

where Xs{x) = u{x)/d{x) and A(x) = u{x)/d{x). At x = .18 we have As(.18) = .454 and 

A(.18) = 1.748 giving rise to A^yi-lS) = —.138 in fair agreement with the experimental 

result -.09 ± .02 ± .025 

The behaviour of A£,y{x) is plotted in Figure 10 together with the experimental point 
measured by NA51 collaboration. 

We consider now the implications for the QPMSR and we begin with the one by Gross 
and Llewellyn-Smith [jl5| 

F3{x) dx = u + d-u-d = 3 . (29) 

Jo 

+.05 

From Table 1 we get for the l.h.s. of Eq.(B^) 2.44 ' in good agreement with the 

— .08 

experimental value 2.50 ± 0.018 (stat.) ± 0.078 (syst.). 
For the l.h.s. of Gottfried sum rule |^ we get 

hu + u-d-d) = .20±. 02 , (30) 

to be compared with .235 ± .026 P]. As long as for the spin sum rules we get 

Am = .62 ±.02 , 

Ad = -.29 ±.04 , (31) 
8 



to be compared with 



Au^ai = 2F = .93 ± .02 , 

Ad^ai = F-D = -.33 ± .02 , (32) 



From Eq.(^) we get 



J„ = -Am + — Arf = .122± .007 , (33) 
1 2 

4 = — An + -Arf= -.030±.010 , (34) 

consistent with the SMC result .136±.011±.011 |[ and the E143 result .129± .004± .009 
g for Ip, and with the E142 resuh -.022 ± .011 [|1|] for /„. 
For the Bjorken sum rule |]13[ we get 



= ^(Au- Ad) = .152±.010 , (35) 

smaller than {l/6)\gA/ Qvl = -209. 

In Table 2. we compare our evaluations of l.h.s. QPMSR with the experiment and 
the prediction of theory without the QCD corrections. 

By comparing the value found for the first momenta (often called by us more prosas- 
tically abundances) of the gas components of the different parton species with the r.h.s.'s 
of Eqs.(|I]) and (0), one finds 

=1.15 ±.01 < M^^; = 1 + F = 1.464 ± .009 , 

=.62 ±.01 <c^L= ^+^~^ = .665±.009 , 

= .53 ± .01 ^ ui^i = l-F = .536 ± .009 , 

dl, =.33 ±.03 ^c^L = ^^^~^ = -335±.009 . (36) 

The different behaviour of with respect to the other valence quark, for which g^^s ~ Ivai, 
may be understood in the framework described here as an effect of Pauli blocking, since 
its levels are almost completely occupied differently from the other valence quarks with 
smaller potentials, as it is also shown by the fact that x('U^) takes the highest value 
allowed. Thus, the interpretation of the defect in Gottfried sum rule as a consequence of 
Pauli principle, disfavouring the most abundant valence parton, u\ seems supported by 
the inequalities (^6]). This interpretation would bring to the very relevant consequence of 



a defect in the Bjorken sum rule. This conclusion is also supported by the good agreement 
with the data of the relationship 

, d 1 , ^ 

" 2 + 2"^ ' ^^^^ 

which implies 

Au = u-d + 2F-l . (38) 

With the abundances found by us Eq.(^) reads 

.53±.01 = .51±.03 . (39) 

A word of caution is welcome for our conclusions on the violation of Bj sum rule, since 
we did not include the effect of QCD corrections in relating the quark parton distributions 
to the structure functions. Also we assumed no polarization for q, being unable to get a 
reliable evaluation of Aq with the present precision for the polarized structure functions 
at small x. Indeed our description of g{{x) and gi{x) is good in terms of Au{x) and 
Ad{x), but our prediction is smaller than the central values of the three lowest x values 
measured by SMC. 

4 Conclusions 

We compared with data the quark-parton distributions given by the sum of Fermi-Dirac 
functions and of a term not contributing to the QPMSR relevant at small x. We obtain 
a fair description for the unpolarized and polarized structure functions of the nucleons as 
well as for the F^lx) precisely measured in (anti) neutrino induced deep-inelastic reactions 
and for q total distribution. The conjectures of previous works on d distributions are 
well confirmed by the values chosen for their thermodynamical potentials. As long as the 
implications for QPMSR the values found for the first momenta of the various parton 
species give l.h.s.'s consistent with experiment. For the fundamental issue of the Bjorken 
sum rule, as advocated in previous works 0, and pS), we get 



Am ^ u-d + 2F -1 , (40) 
Ad > F-D , (41) 

to confirm the suspicion of a violation of Bjorken sum rule related to the defect in the 
Gottfried sum rule. 
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Table 1. 



Parameters 


Previous fit 


Fit BS [17] 


Actual fit 










= 2.47 


A 


.Oo 




+ 09 
2.66 _ 


a 


—.00 


-.646 for ul^i 
-.262 for ul, 


-.203 ± .013 








2.34 +-0^ 
-.06 


A 






0895 +-°^°^ 
■^^^^ -.0084 


OLL 






-1.19 ± .02 








^•^^ -1.59 


X 


.132 


.092 


.235 ± .009 


gas abund. 


x{u^) 


.524 


.510 for ul^i 


1.00 ± .07 


1.15 ± .01 


x{u^) 


.143 


.231 for 


.123 ± .012 


.53± .01 








-.024 


.33± .03 


x{d^) 






200 


.62± .01 


x{u^) 


-.216 




-.886 ± .266 


015 

•^^^ -.009 


x{u^) 


-.141 






// 


= x{d-^) 


// 




- 460 +-°^' 
-.064 


08 

-.02 
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Table 2. 



Sum rule 


Experimental data 


Our fit 


QPM 


GLS 
G 

EJ 1 

[ In 

Bj 


2.50 ± .018 ± .078 \14\ 

.235 ± .026 
.136 ±.011 ±.011 [§ 
.129 ±.004 ±.009 [| 

-.022 ±.011 18] 


2 44 

^■^^ -.07 
.20± .02 

.122 ± .007 
-.030 ± .010 

.152 ± .010 


3 

1/3 
.188 ± .005 

-.021 ±.005 
.209 
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Table captions 



Table 1. Comparison of the values for the parameters of our best fit with the corresponding 



quantities, if any, found in previous analysis H, [117 



Table 2. Comparison of our predictions for the sum rules with the experimental values and 
with the quark parton model (QPM) predictions without QCD corrections. 

Figure captions 

Figure 1. The prediction for -^^(a;) — -F^(x) is plotted and compared with the experimental 
data §. 

Figure 2. The prediction for F2{x) / -^^'(a;) is plotted and compared with the experimental data 
i- 

Figure 3. xgi{x) is plotted and compared with the data 0, 0. 



Figure 4. xg'^{x) is plotted and compared with the data |TH. 



Figure 5. xF-i{x) is plotted and the experimental values are taken from 



Figure 6. xq{x) versus x is shown, the experimental data correspond to [11]. 



Figure 7. The momentum distributions of gas component of q and g's, and of the total liquid 
part are here shown. 

Figure 8. xG{x) versus x is shown, the experimental data correspond to CDHSW [p!9|] . 



SLAC+BCDMS ^ and NMC ^ 



Figure 9. The predicted values for x 
(full line) are compared. 



(dashed line), and l[Fi{x) - F^{x)] 



Figure 10. The asymmetry AoYix) is here plotted, the experimental result is taken from 
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Fig. 4 
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